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Abstract. We summarize recent developments in the study of the origin 
of halo spin profiles and preliminary implications on disk formation. The 
specific angular-momentum distributions within halos in N-body simula- 
tions match a universal shape, M{<j)(xj/{jo+j). It is characterized by 
a power law over most of the mass, and one shape parameter in addition 
to the spin parameter A. The angular momentum tends to be aligned 
throughout the halo and of cylindrical symmetry. Even if angular mo- 
mentum is conserved during baryonic infall, the resultant disk density 
profile is predicted to deviate from exponential, with a denser core and 
an extended tail. A slightly corrected version of the scaling relation due 
to linear tidal-torque theory is used to explain the origin of a typical 
power-law profile in shells, j{M) oc with s ^ 1. While linear theory 
crudely predicts the amplitudes of halo spins, it is not a good predictor 
of their directions. Independently, mergers of halos are found to produce 
a similar profile due to j transfer from the orbit to the product halo via 
dynamical friction and tidal stripping. The halo spin is correlated with 
having a recent major merger, though this correlation is weakened by 
mass loss. These two effects are due to a correlation between the spins of 
neighboring halos and their orbit, leading to prograde mergers. 



1. Introduction 



Disk galaxies that form in cosmological hydro simulations are surprisingly smaller 
and of lower spin than observed disks, posing an "angular-momentum problem" 
(e.g. Navarro & Steinmetz 1997, 2000). A key ingredient in the origin of galactic 
rotation, which can be studied in detail, is the distribution of angular momentum 
in the dark halo within which baryonic infall is assumed to form the disk. We 
describe here several advances in the investigation of the origin of halo angular 
momentum and its distribution, and preliminary implications on disk formation, 
based on high-resolution A^-body simulations and analytic approximations. The 
universal j profile (§2) and the implied disk structure (§3) are presented in Bul- 
lock et al. (2000, BDK+). The issues of tidal-torque theory (§4) are studied by 
Porciani, Dekel, & Hoffman (2000, PDH) and the scaling relation is by Porciani 
&: Dekel (2000, PD). The j profile due to mergers (§5) is by Dekel & Burkert 
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Figure 1. Mass distribution of specific angular momentum and tfie 
corresponding fits by the universal profile for 4 representative halos of 
different /x values. Left: scaled to coincide at the virial radius. Right: 
scaled to coincide at the bending point (jo,Mo). 

(2000, DB). The correlations with merger history (§6) are to be described in 
Wechsler et al. (2000, WBD+). 



2. A Universal Halo Spin Profile 

We use a high-resolution A^-body simulation of the ACDM cosmology (r2m = 0.3) 
inside a 60 /i~^Mpc box, with particles of lO^M^ and force resolution of 2 /i~^kpc 
(ART code, Kravtsov et al. 1997). Spherical halos are identified around density 
maxima using an NFW density profile; they contain a mass M inside a virial 
radius R, defined such that it encompasses a mean overdensity of 340. We focus 
here on 500 isolated halos with M > W^'^h~'^M^ . 

We find that the mass distribution of specific angular momentum j is well 
fit by a universal profile: 

M{<j) = fiMj/{jo+j), /i>l. (1) 

This is roughly a power law for j < jo, which flattens off if j > Jq. The shape 
parameter /i correlates with the fraction of mass in the power-law regime, which 
always contains more than half the mass. 

One can replace jo with a robust spin parameter A' = J/{y/2MVcR), where 
= GM/ R. It coincides with the standard spin parameter A for an isothermal 
sphere and for an NFW profile with concentration C ~ 10. The pairs (/i,jo) 
and (//, A') are related via 

jo 6(/i) = V2V,R A', b{^l) = ln(l - ^-1) - 1 . (2) 

Either pair can fully characterize the j profile. Fig. 1 shows a few examples. 
Fig. 2 (left) shows a certain correlation between ^ and A', with a correlation co- 
efficient r = 0.23, meaning that high-spin halos tend to have power-law profiles. 
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Figure 2. Left: Correlation between angular-momentum parameters 
/i and A'. The symbol size is inversely proportional to the error. Re- 
gression lines are shown; the correlation coefficient is r = 0.23. Right: 
Disk surface-density profile as implied by the halo M{<j) with j con- 
servation. Shown for comparison are an exponential disk and the disk 
resulting from a uniform, solid-body rotating halo. Note the density 
excess at small and large radii. 



We find that the direction of the angular-momentum vector tends to be 
aligned throughout the halo; for more than half the mass, j and jz typically differ 
by less than a factor of 2. The spatial symmetry of j tends to be more cylindrical 
than spherical; in 80% of the halos, the mean value of j (over longitude) decreases 
by less than a factor of 2 when moving away from the equatorial plane along 
lines parallel to the spin axis. 



3. Disk Structure: Core and Tail 

A simple model can convert the halo M(< j) into a disk density profile. Assume 
that the gas initially traces the j distribution of the halo, and that j is conserved 
during the gas contraction into a rotation-supported disk of mass fM. Then, 
using Eq. 1, the disk mass inside radius r is oc j(r)/[jo + j{r)\. Under the 
assumption of an isothermal halo with a flat rotation curve, one has j(r) = rV^, 
and then the disk surface-density profile is 

_ //iM rd _ V2\'R 

27r r(rd + ' '''^ ' ' ^ ' 

The characteristic length scale is based on Eq. 2. 

A more accurate numerical calculation using an NEW profile, that self- 
consistently solves for the disk circular velocity at r and takes into account the 
adiabatic contraction of the halo in response to the baryonic infall, leads to the 
typical profile shown in Eig. 2 (right). Compared to an exponential disk (which 
is similar to the profile obtained from a uniform, solid-body rotating halo), our 
resultant profile has a core and tail excess, with somewhat higher densities both 
at small and large radii. Since observed (stellar) disks are closer to exponential. 
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this result may either indicate that j is not exactly preserved during the collapse, 
or it may require some additional viscous processes (e.g. Lin & Pringle 1987). 
In some cases, the obtained slow-rotating core may be associated with a bulge. 

The extended disk, on the other hand, may be useful in explaining the damped 
Lyman-alpha systems observed at high redshift (Mailer et al. 2000). 

Our result for the limiting case of j conservation clearly sharpens the big 
"angular-momentum problem" , because j transport from gas to dark matter and 
from inside out can only worsen any agreement with observed disks. It seems 
that one needs to appeal to some process of gas removal from subclumps before 
they merge into the halo center, and thus reduce j transport. 

4. Spin Profile based on Tidal Torque Theory 

We find in the simulations that the halo j profile in spherical shells is roughly 

a power law, j{M) oc M®, with s distributed like a Gaussian: s = 1.3 it 0.3 (see 
an early indication by Barnes &: Efstathiou 1987). Here j is the specific angular 
momentum in a shell of radius r, and M is the mass in the sphere encompassed 
by this shell. This profile can be related to the distribution M(< j) under the 
assumptions of spherical symmetry for the density and cylindrical symmetry for 
j. Linear tidal-torque theory (TT) provides one hint for the origin of this power 
law, as follows. 

TT (e.g. White 1984) implies that the j gained by a halo at time t before 
its turn-around is 

Ji{t) = a{tfD{t)eijkTjiIik, (4) 

where the time growth is from some fiducial initial time t[, Iik is the inertia tensor 
of the proto-halo at ti, and Tji is the tidal tensor at the halo center, smoothed on 
the halo scale. This is based on assuming the Zel'dovich approximation for the 
velocities inside the proto-halo, and a 2°^^-ordcr Taylor expansion of the potential 
(see PDH). We show elsewhere (PD) that the standard scaling relation of TT 
should be slightly modified; it should read 

joctcfT(M)M2/3 , (5) 

where is the turn-around time of the halo (assuming that the relevant cos- 
mology at tc is still EdS) and a{M) is the rms density fluctuation on scale M 
at tj. 

To obtain the spin profile we apply Eq. 5 shell by shell. The turn-around 
time of each shell is determined by 5{M)D{tc) ~ 1, where 5{M) is the mean 
density inside M. The densities within the different shells are correlated, such 
that, for a Gaussian field, the typical density fluctuation profile about a random 
point scales like 5{r) oc ^(r), where ^(r) is the linear two-point correlation func- 
tion (Dekel 1981). This is accurate to a few percent also around a high peak 
(Bardeen et al. 1986, Fig. 8). Thus, for a power-law power spectrum P/^, oc A;" 
and a flat universe, we obtain for the j profile within each halo 

i(M) oc M2/3+(3+n)/3_ (g) 

This implies 1 < s < 4/3 for — 2 < n < — 1, the range appropriate for large 
galactic halos in a CDM spectrum, in pleasant agreement with our finding for 
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Figure 3. Left: Angular-momentum profile in spherical shells due 
to a minor merger of a satellite into an NFW halo, assuming circular 
orbits and 3 different recipes for tidal mass stripping. The dotted lines 
are of slopes 1 and 2 for reference. Right: Halo spin correlations in a 
cosmological A^-body simulation. Shown are the averages of cos(angle) 
between the spins and between spin and orbit for all pairs of halos 
separated by r. 

the simulated halos. A similar analysis, using the accretion rate into halos 
based on the Extended Press-Schechter approximation, also explains a weak 
anti-correlation detected in the simulations between s and halo mass (BDK-I-). 

Using simulated halos (in GIF simulations), PDH find that TT predicts an 
amplitude for the halo spin in the right ball-park. This confirms the relevance 
of TT to the j distribution in halos. On the other hand, PDH find that the 
predicted spin direction is only weakly correlated with the true direction, with 
an average deviation of 53°. This refiects large inaccuracies in the smoothing of 
the tidal tensor and the truncation of the Taylor expansion of the potential. This 
relative failure of TT in predicting the spin direction should be born in mind 
when trying to use TT for the analysis of spin-spin correlations and weak-lensing 
phenomena. 

5. Spin Profile due to Mergers 

A simple toy calculation for how j is being deposited in the halo during a merger 
provides another hint for the origin of the j profile. Consider a fixed halo and 
an incoming satellite of mass profiles M(r) and m{i) respectively. Assume that 
as it moves in, the satellite is losing mass outside a tidal radius it, determined 
at r by 

m{it)/it = 4/u(r)/r^ , /i(r) = [2M(r) -rdM/dr] . (7) 

For isothermal halos this implies m(r) oc M{r) oc r. Now assume that the 
satellite is spiraling in due to dynamical friction in circular orbits, and that the 
mass and j are deposited locally. Then j{r) is obtained by averaging over shells, 

47rrV(r) j(r) = m(r) ^ ^ + ^ rV^ir) . (8) 
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The j transfer is due to both the slowdown of the satelhte by dynamical friction 
and the direct tidal stripping of mass. For an isothermal halo, we obtain 

j{M) cx M . (9) 

Fig. 3 (left) shows the j profile for an NFW halo swallowing a satellite that 
is loosing mass either according to the recipe m(r) oc M(r) or based on Eq. 7 
assuming that the satellite has an isothermal profile, compared to the case of 
no mass loss. The general shape of j{M) in the outer decade of mass, when 
realistic mass loss is considered, is reminiscent of the profiles detected in the 
cosmological simulation. 

Full iV-body simulations of mergers, spanning a range of halo and collision 
parameters, confirm the robust production of such power-law like profiles (DB). 

The resulting j{M) from a sequence of minor mergers would be a sum of 
similar contributions, each projected onto the direction of the total net angular 
momentum (perhaps determined by the most major merger). 

Since halo accretion histories typically reflect some combination of relatively 
quiescent mass accretion and more pronounced mergers, the j profiles of indi- 
vidual halos may reflect a combination of the two processes explored above. We 
see that the profiles produced by each of the two processes are, in general, of 
similar shape to the profiles measured in the simulations. 

6. Spin and Merger History 

Further investigations of halos in the ART and GIF simulations (WBD-I-) reveal 

interesting correlations between halo spin and merger history. For example, halos 
that had a major merger in the last Gyr tend to have higher spins, (A) = 0.055, 
in comparison with the average of 0.035 for all halos. This is hardly surprising 
because most of the orbits of merging halos are expected to be of large impact 
parameter and therefore to bring in much angular momentum. 

On the other hand, we find that major mergers also result in effective loss 
of j, typically by a factor of 2 of the original orbital j. This spin loss is 
associated with mass loss of the weakly bound outer layers of the halos which 
are preferentially j-rich, Eq. 1. This spin loss somewhat weakens the tendency 
for high spin in merger products. 

These effects may be correlated with the fact that major mergers tend to 
be prograde, with the spins of the two merging halos oriented roughly parallel 
to each other and to the angular momentum of their mutual orbit. Fig. 3 (right) 
shows the average cos{6) between the two spins and between the spin and orbit 
for pairs of halos at a given separation r. There is a significant alignment at 
r < 400^~^kpc, growing as a power law towards smaller distances. Prograde 
encounters indeed explain both the high spin and the high spin loss. Being 
prograde they are more effective in transferring energy from the orbit to internal 
motions and thus making the typical impact parameter for mergers even larger. 
This, plus the coherency of progenitor spins and orbital j, results in a high spin. 
But at the same time a prograde encounter causes effective stripping of j-rich 
mass. 

High spin due to major mergers seems not to go naturally with the common 
wisdom that such mergers produce low-spin bulges and elliptical galaxies. How- 
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ever, it is the halo properties at z ~ 1 (say) rather than today that are hkely to 
determine the j content of the gas which later falls in to form the disk. Then, 
it is the halo merger history prior to 2: ~ 1 that is more relevant for the galactic 
spin, while later mergers govern the bulge-to-disk ratio. The high-spin halos at 
z ~ 1, those that were assembled and virialized not long before z ~ 1, may be 
the hosts of high-spin disks that are observed as low-surface-brightness galaxies 
(Dekel et al. 2001). 
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